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Abstract 

For ordinary hcrmitian Hamiltonians, the states show the Kramers degeneracy when the 
system has a half-odd- integer spin and the time reversal operator obeys O 2 = —1, but no such 
a degeneracy exists when 2 = +1. Here we point out that for non-hermitian systems, there 
exists a degeneracy similar to Kramers even when 2 = +1. It is found that the new degeneracy 
follows from the mathematical structure of split-quaternion, instead of quaternion from which 
the Kramers degeneracy follows in the usual hcrmitian cases. Furthermore, we also show that 
particle/hole symmetry gives rise to a pair of states with opposite energies on the basis of the 
split-quaternion in a class of non-hcrmitian Hamiltonians. As concrete examples, we examine 
in detail N x N Hamiltonians with N = 2 and 4 which are non-hermitian generalizations of 
spin 1/2 Hamiltonian and quadrupole Hamiltonian of spin 3/2, respectively. 
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1 Introduction 

The original observation between time-reversal (TR) invariance and statistical mechanics is 
traced back to Dyson who pointed out that TR operator is naturally incorporated in the 
algebra of quaternions if the system has a half-odd-integer spin and Q 2 = —1 [I]. He showed that 
the Kramers degeneracy comes from mathematical structures of quaternion, and its statistical 
properties are described by the symplectic group. On the basis of these, Avron et al. explored 
topological properties of fermionic systems with TR symmetry [2], and the second Chern number 
was introduced as an extension of the TKNN topological number of quantum Hall effect (3J 0] . The 
TR symmetry and the resultant Kramers degeneracy also play a central role in recent developments 
of the quantum spin Hall effects [6] and topological insulators [U [9]. Indeed, the Kramers 
degeneracy enables us to introduce a new class of topological numbers characterizing these phases. 
The mathematical structures of the topological insulators have been studied in Refs. [10l 112]. 

Meanwhile, if the system has an integer spin and G 2 = +1 such as boson systems and systems 
with even number of electrons, we have no such a Kramers degeneracy. Correspondingly, its 
topological structure is rather simple and the Hamiltonian has a real structure in general. However, 
such a consequence changes if we allow non-hermiticity of Hamiltonians. Indeed, as is shown 
below, there is generally a degeneracy similar to the Kramers even when 2 = +1 in a class of 
non-hermitian Hamiltonians. 
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Although we usually suppose hermiticity of Hamiltonian, non-hermitian Hamiltonians also 
have applications to interesting problems such as open chaotic scattering [13], dissipative quan- 
tum maps [13], and derealization of pinned vortices in superconductors [T5]. We also have 
non-hermitian Hamiltonians as effective theories of hermitian systems. Moreover, non-hermitian 
Hamiltonians might be meaningful themselves if a kind of TR symmetry such as VT symmetry 
|16j or pseudo-hermiticity [T7] is imposed. They are a part of the motivations that we pursue the 
present work. 

We investigate TR symmetry with 2 = +1 in non-hermitian Hamiltonians. From a gen- 
eral argument, it is shown that such symmetry is naturally incorporated in the algebra of split- 
quaternion, instead of quaternion. (See also related work [18^ [T9] in VT symmetric quantum 
mechanics.) Then a new kind of degeneracy is obtained from structures of split-quaternion. As 
concrete examples, we examine N x N non-hermitian Hamiltonians up to A = 4. The structure 
of split-quaternion is identified in these Hamiltonians, and we find that it has a close similarity 
to the quaternion structure of the spin 1/2 Hamiltonian and quadrupole Hamiltonian of spin 3/2. 
Furthermore, it is shown that the particle/hole symmetry also gives rise to a pair of states with 
opposite energies (E, —E) in a class of non-hermitian Hamiltonians. Random matrix classification 
of the non-hermitian models is also provided. 

The organization of the paper is as follows. In Secj2j a generalized Kramers degeneracy 
in pseudo-Hermitian quantum mechanics is discussed. We point out relations between split- 
quaternion and TR operation for G 2 = +1, and show the existence of generalized Kramers degen- 
eracy in pseudo-hermitian systems. We also argue relations between the particle/hole symmetry 
and split-quaternions. The split-quaternion structure of integer spin systems is clarified, too. In 
SecJU we show how the generalized Kramers theorem in pseudo-hermitian systems is incorporated 
in the non-hermitian random matrix classification. As a concrete example of pseudo-hermitian 
model with particle/hole symmetry, SU (1, 1) model is introduced in SecJU and basic properties of 
the model are investigated. In SecJSl we argue properties of the SO(3, 2) model with time reversal 
symmetry 2 = +1 as a simple exemplification of the generalized Kramers degeneracy. It is also 
shown that the SO(3, 2) model is realized as a SU(1, 1) quadrupole model with SU(1, 1) spin 3/2. 
Secj6]is devoted to summary and discussions. 

2 Generalized Kramers degeneracy and split-quaternion 

2.1 Split-quaternion and time-reversal symmetry 

Let us start with a brief review of the split-quaternion. The split-quaternion |20] is a variant 
of the quaternion [21 j which is written as 



with real numbers (w,x,y,z) in the basis (l,i,j,k). The algebra of the basis for the split- 
quaternion is different from that for the quaternion, and it is given by 



q = w + xi + yj + zk, 



(1) 



ij = k= -ji, jk 



kj, hi = j 



ik, i 2 = -1, j 2 = 1, k 2 = 1. 



(2) 
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Note that j 2 = 1 and k 2 = 1, not j 2 = — 1 and k 2 = — 1 as in the quaternion case. 

It has been known that the structure of the quaternion naturally arises in time-reversal (TR) 
invariant systems. The TR operator is antiunitary, and anticommutes with i: 

= UK (3) 

where U is a unitary operator and K complex conjugates everything to its right. For systems with 
an integer spin, we have 2 = +1, while for systems with a half integer spin, 2 = —10 In the 
latter case, the TR invariance results in the structure of quaternion for the Hamiltonian. Then, 
what is a natural mathematical structure in the former ? 

The answer is the split-quaternion. The TR operator is antiunitary 

Qi = (4) 

By identifying j and k with and iO, respectively, one finds a correspondence between the triplet 
of the TR algebra and the split-quaternion, 

(i,Q,iQ) o (i,j,k). (5) 

Thus, the split-quaternion also fits into the TR symmetry with 2 = +1. 

In spite of the argument above, it has been known that there is no such a split-quaternion 
structure in the TR invariant Hamiltonians with 2 = +1. The Hamiltonian supports only a 
real structure instead [I]. This is because usually the Hamiltonians are supposed to be hermitian. 
This implicit assumption makes the split-quaternion into a real number. Nevertheless, physical 
phenomena are not always described by hermitian Hamiltonians. Non-hermitian Hamiltonians also 
have interesting physical applications. Then, if we consider a class of non-hermitian Hamiltonians, 
the hidden split-quaternion structure becomes evident, as will be shown in the following sections. 

2.2 Pseudo-hermiticity 

A non-hermitian Hamiltonian H is called pseudo-hermitian j!7j . when it satisfies pseudo- 
hermiticity 

H^ = r]Hr]-\ (6) 

where r\ is a hermitian operator referred to as the metric operator. For example, consider a 
nonunitary transformation G on a hermitian Hamiltonian Hq, then H = GHqG^ 1 is not hermitian, 
but pseudo-hermitian, 

ift = GTt-ifljjGTt = ($-icr x HG&, (7) 

with a metric operator 7] = (GG^)^ 1 . 

1 The action of the TR operator is two, i. e. the system comes back to the original if we apply Q twice. Thus Q 2 
should be a phase e ia , which implies U = e ia U T and U T = Ue ia . This yields U = Ue 2ia , so the phase is e ia = ±1. 
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The reason why rj is called the metric operator is that the time-independent inner product of 
a state is given by a metric rj. For a non-hermitian Hamiltonian, there generally exists a set of 
states, \4>) and \<p)) that satisfy [22] 

i^) = H\0), ij^p»=jft|p». (8 ) 

The time independent inner product is constructed as {(j)\<p)), as shown by 

For a pseudo-hermitian Hamiltonian ff, © leads to 

i^) = rfT,\<l>), i|r,-V)) = HrrV)). (10) 

Thus, we also have additional time-independent inner products, (cf)\r)\(f)) , (((/3|ry _1 |(/5)). In the fol- 
lowing, we mainly use {<fi\(p)) as the inner product unless explicitly written. 

In order for the pseudo-hermiticity to be consistent with the TR symmetry, the condition ([6]) 
should be commutative with the TR operation. This leads to r/* = ±U^r]U , that is Qrj = ±7/0. 
Therefore, possible metric operators are classified into two: The first one satisfies [rj, G] = 0, 
and the second {rj, G} = 0. Here we note that the latter case is proper for only non-hermitian 
Hamiltonians. For hermitian Hamiltonians, we have rj = 1. Thus r/ commutes with G rather 
trivially. 

2.3 Generalized Kramers degeneracy 

Let us assume two conditions: one is the TR symmetry with G 2 = +1 

[H,e] = o, (ii) 

and the other is the anticommutation relation 

{ V ,Q} = 0. (12) 

Let \4> n ) be eigenstates of H, 

H\cf> n ) = E n \<j> n ). (13) 
Then the corresponding eigenstates \<ft n )) of H\ 

H^ n )) = E* n \<p n )), (14) 

which satisfy 

t>n) = 5nm- (15) 



■"n\Ym — \\™ 
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The eigenstates \4> n ) and \<j> n }} satisfying (|15p are known as the bi-orthonormal basis [23|, I22j. By 
using the pseudo-hermiticity ©, is rewritten as 

Hrj~ 1 \<p n )) = E* n r ] - 1 \<p n )). (16) 

We apply from the left to both sides of (|16p to have 

He V - 1 \<p n }) = E n er ] - 1 \ ( / )n }}, (17) 

where on the left-hand side, we utilized the time reversal invariance of the Hamiltonian (|lip . Thus 
\(f) n ) and ®rj~ l \<j) n )) have the same eigenvalue E n . Therefore, if they are linearly independent, we 
have degeneracy in the eigenvalues of H. Note 

((</g07rV4 = ((eVVniec/a = i^- 1 ®^)) = -((0 n |e?rV4- (is) 

In the first equation, we have used the antiunitary property of 0, and the second equation follows 
from 2 = +1 and the hermiticity of rf 1 . In the third equation, the anticommutation relation 
between G and rj~ l was utilized El- Thus we have ((4> n \®ri~ 1 \(j) n }} = 0. On the other hand, 
((4>n\4>n) = 1 from ([15]) . Therefore \4> n ) and ®r]~ 1 \(f) n )) are linearly independent [§. As a result, we 
have two-fold degeneracy in eigenstates of H, which is the generalized Kramers degeneracy. 

In general, the generalized Kramers partner 07? _1 |</> n }) is not coincident with 0| </>„,). Actually, 
unlike the TR symmetry with 2 = —1, TR symmetry with 2 = +1 does not imply that \4> n ) 
and @\4> n ) are linearly independent. Nevertheless, we can say that if eigenvalues of H are real, 
the generalized Kramers partner is essentially the same as ®\4> n ). 

To see this, let us consider the eigenstate \<f> n ) satisfying (fl3j) . Then the pseudo-hermiticity 
leads to 

H^\4> n ) = En^n). (19) 
Therefore, r]\4> n ) can be expanded as 

r]\4>n) = \4>m))Cmn, (20) 
m 

where the sum is taken for |^ m ))'s satisfying (I14p and (|15p with E^ = E n . (Note that if there is 
a degeneracy in the spectrum, we may have multiple such |(/> m ))'s.) Applying (4> m \ from the left, 
we obtain 

Cmn = {<j>m\v\4>n) ■ (21) 

Because of the hermiticity of r/, c mn is hermitian for the indices m and n. Thus it can be 
diagonalized by a unitary matrix G 

y^C^QfcCfcn = A m <5 mra , (22) 
Ik 

2 For usual hermitian Hamiltonians, we have 77 = 1 in general. Thus, the anticommutativity (|12jl does not hold. 
This is the reason why even if a hermitian Hamiltonian has TR symmetry with O 2 = +1, there is no generalized 
Kramers pair. 

3 Though \4> n ) and ©77" 1 are linearly independent, they are not orthogonal, i.e. (0„|O7; _:L \4> n }) 7^ 0, in 
general. 
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with real A m . The eigenvalue A n is not zero because c mn is invertible. Thus taking the following 
new bi-orthonormal basis 

l^}-El^> G »/^ l« = EW G ™VW. (0nlO=<W, (23) 
m m 

we have 

7 ? |<} = S gn(A n )|<)). (24) 

Now suppose that E n is real, then E m for |(/> m )) in (|20|) is also real and coincides with This 
yields that all |</> m )'s in the right hand side of the first equation in (|23p have the same energy E n . 
In other words, \c/J n ) remains to be an eigenstate of H with the eigenvalue of E n . Applying O77 -1 
from the left to both sides of (|24p . we find that the (generalized) Kramers partner @i]~ 1 \(p' n }} is 
the same as @\4>' n ) up to an irrelevant overall sign, 

eK) = S gn(A n )e ? rVn>>- (25) 



2.4 Particle/hole symmetry and split-quaternion 

In addition to the TR invariance, we can have the particle/hole symmetry C which is antiuni- 
tary. Here we briefly see the split-quaternion structure of particle/hole symmetric system. 
We say that a Hamiltonian H has the particle/hole symmetry C if it satisfies 

CHC' 1 = -H. (26) 

If we write C as C = TK with a unitary operator T, (j26|) is recast into 

YHY ] = -H*. (27) 

One can show that C 2 = ±1. 

In a manner similar to the TR symmetry, we have the split-quaternion structure if C 2 = +1. 
The correspondence between the particle/hole symmetry and the split-quaternion is 

(i,C,iC)^(i,j,k). (28) 

When H is pseudo-hermitian, = i]Hr]~ l with {C, rf\ = 0, we find that eigenstates of H are 
paired with eigenvalues (E n , —E n ). Consider 

H\(j) n ) = E n \(t> n ), 

fft|^»=KI&,», (29) 

with ((4> n \4> m ) = {4> m \4>n)) = $mn- It is found that \<f> n ) and Ct? -1 !^™)) have the eigenenergies E n 
and —E n , respectively. Then, if -q satisfies {rj,C} = 0, we can show that \<p n ) and Cr/ _1 |^ n )) are 
linearly independent for any E n , in a manner similar to Sec l2,31 Thus, the eigenstates of H are 
paired. In particular, if we have a zero energy state with E = 0, then it should be degenerated. 
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Note that the particle/hole symmetry itself implies that if \<j) n ) is an eigenstate of H with 
eigenenergy E n , then C\4> n ) is the one with —E*. For a non-hermitian Hamiltonian, however, this 
does not always mean additional pair of states. Indeed, if E n is pure imaginary, then E n is the 
same as — -E 1 *, and \<p n ) and C\(p n ) can be the same. We can also say that if E n is real, C\(p n ) is 
essentially the same as C^ _1 |(/) n )): When the eigenenergies E n are real, we can take the basis ([231) 
as the eigenstates of H, which leads to 

C|^) = sgn(A„)C^V / n))- (30) 

Thus C\4>' n ) and Cry -1 !^)) coincide with each other up to a sign factor. 

Formally we can treat the particle/hole symmetry as the TR symmetry by redefining H — > iH. 
In this case, however, the pseudo-hermiticity is replaced by "pseudo-anti-hermiticity" , 

H ] = -riHrj- 1 . (31) 



2.5 Example: 2x2 matrix 

In this subsection, we will see the split-quaternion structure in a concrete example. Consider 
a 2 x 2 matrix, as the simplest nontrivial Hamiltonian. In general, by using the 2x2 unit matrix 
I2 and the Pauli matrices <Ji (i = x, y, z), any 2x2 matrix can be written as 

H = hl 2 + Yl hi(Ti > ( 32 ) 

i=x,y,z 

with complex numbers h and h l (i = x,y,z). Then suppose that H is invariant under the TR 
symmetry = UK with 2 = +1. 

[H, 6] = 0. (33) 

2 = +1 implies that U is a symmetric (unitary) matrix, U = U T . Following Ref.pQ, U can be 
U = I2 in a proper basis of the Hamiltonian. The TR invariance yields 

h* = h, h x * = h x , h y * = -h y , h z * = h z , (34) 

thus, we obtain 

H = wl 2 + xa x + yiay + za z , (35) 

with real numbers, w, x, y, z. The split-quaternion structure of H is evident if we notice the 
following identification between the Pauli matrices and the basis for the split-quaternion, 

(ia y ,a x ,a z ) -H- (i,j,k), (36) 

which reproduces the algebra ([2]). Thus the TR invariant Hamiltonian ()35[) is a split-quaternion. 

Let us now impose the pseudo-hermiticity. To satisfy {0, rf\ = 0, the hermitian matrix 7/ 
should be pure imaginary, rf = —77. Thus it can be written as ij = ca y with a real number c. 
If H is pseudo-hermitian ([6]), we obtain x = y = z = 0. Therefore, the Hamiltonian becomes 
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H = wl%. While this Hamiltonian is rather trivial, we have a two- fold degeneracy. This can be 
considered as the generalized Kramers degeneracy explained in Sec l2.31 In this case, any column 
vectors \<fi) = (a,b) T are eigenstates of H. The corresponding \(f>)) satisfying {4>\4>)) = 1 is 



1 / a 



Ir// |a| 2 + |6| 2 \b )' 
Thus the generalized Kramers partner, Qrj" 1 ^)), is given by 



(37) 



e ""^» = ~WTW" y ( t- ) - WTW ( ) ' (38) 

One can easily check that \4>) and Qr]~ l \cf>)) are linearly independent if \a\ 2 + \b\ 2 ^ 0. 
Next, impose the pseudo-anti-hermiticity which implies w = in ([35]) . so 



H = xa x + yicjy + za z . (39) 
The eigenvalues are E± = i-^/x 2 + z 2 — y 2 . The corresponding eigenstates \4>±) are given by 

I0±) = c± ^ ± x ^ + _^ , s = V^ + ^-y 2 , (40) 

where c± are constants. In accordance with the general argument in Sec l2.41 the eigenstate with 
the eigenvalue E is paired with the one with — E. 

We can also check that if we suppose the hermiticity of H, instead of the pseudo-hermiticity 
or pseudo-anti-hermiticity, the split-quaternion structure is replaced by the real structure, y in 
being zero, thus H reduces to a 2 x 2 real symmetric matrix. 



2.6 Integer spin systems 

In this subsection, we will see the split-quaternion structure of integer spin systems in which 
2 = +1. The spin, Si (i = x, y, z), changes its sign under the TR transformation 

es^e- 1 = -si. (4i) 

Write G = UK, and the condition (|4ip is written as 

US X U- X = -S x , USyU- x =S y , US z U~ l = -S z , (42) 

where we have assumed the standard matrix realization of Si in which only S y is complex and 
pure imaginary. Then, U is given by@ 

U = e iwS \ (43) 

and 

9 2 = e 2iwS y = {-l) 2S , (44) 



4 The form of TR operator depends on the physical meaning of the operator. For instance, when Si denote 
"isospin" labeling two different energy levels, the TR operator is simply given by = K . 
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where 5 represents the magnitude of the spin. Thus, for integer 5, 2 = +1, while for half-integer 
5, 2 = — 1. For low spins, 



5 = 1/2 : G = ia y K, 

( 1 

9 = 



S = 1 



5 = 3/2 : G 



V 



-1 | K, 

1 



lOy 

ia y 



K. 



(45a) 
(45b) 

(45c) 



Let us first consider the 5 = 1 system. The TR operator is given by (|45b|) . Since U is real 
symmetric matrix, it can be diagonalized by an orthogonal matrix O, 

( 1 
0-10 
0-1 



6 = 



(46) 



V 



Then is written as 



= OVV t O t K = (OV)K(OV) 



(47) 



with V = diag(l, i, i). Therefore, performing the unitary transformation OV on the basis, is 
recast into = K. 

Let us now impose the TR invariance on the Hamiltonian. In general, the Hamiltonian H 
is written as 

h + h x a x + ih y o y + h z a z a T \ 

b c r 



H 



(48) 



where h and h l (i = x, y, z) are complex numbers, a and b are two component complex vectors, 
and c is a complex number. Take [H, 0] = with = K, then all the elements of H are real. In 
particular, from the correspondence (|36p . this implies that the left upper part of H is given by a 
split-quaternion, h r = w + xa x + iya y + za z with real coefficients w, x, y, z. 
Next, consider the 5 = 2 case, in which S y is given by 



and the corresponding TR operator becomes 







(o 


-2 













2 





-V6 













v/6 





-V6 













VQ 





-2 




\o 








2 





/ 



(49) 



(o 











1 \ 











-1 











1 











-1 











v 1 











o ) 



K. 



(50) 
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In a manner similar to the above, we have = K by choosing the basis properly. The TR 
invariant Hamiltonian is given by 



/ hf 2 aj 

H = h,2i 0-2 

\ bi b 2 c 



(51) 



where hfj (I, J = 1,2) are split-quaternions of the form 



I j = wij + xij<j x + iyijUy + z u a z , 



(52) 



with real coefficients wij,xu,yu,zjj. Here a/ and bj (I = 1,2) are real two row vectors, and c 
is a real constant. 

In general, for S = M with integer M, we can always choose the basis in which is given by 
G = K. Then a TR invariant Hamiltonian H is given by 



where hfj (I, J = 1, • • • ,M) are split-quaternions, and aj and b[ (I = 1, • • • ,M) are real two 
row vectors, and c is a real constant. 

Note that if we impose the hermiticity on the Hamiltonian (j53[) , H reduces to a real symmetric 
Hamiltonian. Thus, our result here is consistent with the known result that the hermitian TR 
invariant systems with integer spins belong to the orthogonal ensembles [lj [23] . 

3 Random matrix classification 

The idea of random matrix ensembles, pioneered by Wigner and Dyson, is based on classifying 
classes of matrices by discrete symmetries (e.g. see |24j). Altland and Zirnbauer established the 
hermitian random matrix theory in the context of superconductivity \lf>\ I26|. which contains 10 
classes. The random matrix classification was also applied to topological insulators |llj . 

Bernard and LeClair extended the random matrix classification for non-hermitian matrices 
|27j . In this section, we apply the arguments in Secl2]to the framework of the random matrix 
classification. 

3.1 Non-hermitian random matrix classification 

Following Rcf. [27], let us consider discrete symmetries on non-hermitian random matrices. 
Suppose that the discrete symmetries are implemented by unitary transformations, and the system 
comes back to the original up to a phase factor if they are applied twice. Then there are four 




(53) 
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possible transformations on a non-hermitian random matrix H d: 



K 


sym. : 


H = kH*k~\ kk* = ±1 


Q 


sym. : 


H = e q qHU~\ gV 1 = 


c 


sym. : 


H = e c cH T c~ 1 , c T c~ l = 


p 


sym. : 


H = -pHp' 1 , p 2 = 1, 


i.e. 


e c = ±l, 


and e q = ±1, and k, q, c, < 



(54a) 

1, (54b) 
±1, (54c) 

(54d) 



kk* = 1, qq^ = 1, cc' = 1, pp^ = 1. (55) 

Demanding that the transformations (|54ap - (|54dp commute, we have 

q* = ±k~ 1 qk^\ k T c^kc* = ±1, p* = ±k~ 1 pk, 

q T = ±c^g _1 c, q = ±pqp\ c = ±pcp l . (56) 

We refer to K symmetry as the TR symmetry. We can add the minus sign to the right hand 
side of the first equation in (|54ap by redefining H —> iH. Thus we can also consider K symmetry 
as the particle/hole symmetry. 

Q symmetry corresponds to the pseudo-hermiticity (e q = 1) or the pseudo-anti-hermiticity 
(e q = —1), defined in the previous section, by identifying q with r/ _1 . We note that the corre- 
spondence is not one-to-one. While q is a unitary operator, r] is not always. (Both r/ and q are 
hermitian.) Thus Q symmetry is a part of the pseudo-(anti-)hermiticity. 

In the case of hermitian matrices, K symmetry is nothing but C symmetry. Thus, one often 
refers to C symmetry as the TR symmetry if e c = 1 or particle/hole symmetry if e c = —1. For 
non-hermitian matrices, however, they are different. Thus, in this paper, we do not refer to C 
symmetry as the TR symmetry or particle/hole symmetry. C symmetry is obtained by combining 
K and Q symmetries. Specifically, from K (|54ap and Q (|54bp symmetries, C symmetry (|54cp is 
obtained with c = kq* (up to a phase factor) and e c = e q . 

Finally, P symmetry is called the chiral symmetry in literatures. 



3.2 Random matrix classification and split-quaternion 

Write 

k = U, (57) 
then (I54ap gives the TR symmetry with = UK as 

[e,H] = o, e 2 = ±i. (58) 

Next write 

q = if 1 , (59) 

5 The terminology "symmetry" is usually used for the operations that commute with Hamiltonian, but in this 
section, "symmetry" refers to the operations (|54a[) - ()54dl) . 
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then (|54bp reads 

H ] = e^Hrf 1 , ?? f = 77. (60) 

Thus Q symmetry with e q = 1 as the pseudo-hermiticity, and with e q = — 1 as the pseudo-anti- 
hermiticity. 

If the first equation of (|56|) holds, a system has both K and Q symmetries. In terms of and 
77, the commutativity between K and Q, i.e. q* = ±k~ 1 qk^~ l , is written as 

G? ? =F 776 = 0. (61) 

Thus K symmetry and Q symmetry are equipped with all the properties used in the arguments 
in the previous section. 

The arguments in the previous section lead to the following. 

1. When a non-hermitian matrix has K symmetry with kk* = 1, the matrix supports the 
split-quaternion structure. 

2. If the non-hermitian matrix also has Q symmetry with e q = 1 and q* = —k~ 1 qk^~ l , at the 
same time, each eigenvalue of the non-hermitian matrix has two-fold degeneracy. 

3. If the sign of Q symmetry is minus, i.e. e q = — 1 (and if q and k satisfy q* = —k~ 1 qk<~ 1 ), 
each eigenstate with the eigenvalue E of the non-hermitian matrix has a partner state with 
—E. 

The second result is nothing but the generalized Kramers degeneracy in Sec J2.3( and the last one 
comes from particle/hole symmetry arguments in Sec J2.4l 

4 Random matrix class of the 577(1,1) model 

As a concrete realization of the statements in Sec l3.2t we introduce the SU (1, 1) models. First 
consider K symmetry realized by k = a x with kk* = 1. Although k can be k = I2 by taking a 
proper basispQ as mentioned in Sec. 12.51 and explicitly shown below, the present form of k = a x 
is convenient to see the SU(1, 1) structure. 

Any arbitrary 2x2 matrix can be expanded by 2x2 unit and the Pauli matrices 

H = hl 2 + hi(J i> ( 62 ) 

i=x,y,z 

where h and h l (i = x, y, z) stand for complex parameters. The imposition of K symmetry specifies 
h and h l as 

h* = h, h x * = h x , h y * = h y , h z * = -h z . (63) 

Thus, h, h x , h y are real parameters, while h z a pure imaginary parameter. Thus, the Hamiltonian 
is rewritten as 

H = wl 2 + xa x + ya y + iza z , (64) 

with real parameters w, x, y, z. We further impose Q symmetry with q = a z . There are two types 
of Q symmetry, corresponding to e q = ±1. 
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4.1 e q = +1 : pseudo-hermiticity 

In this case, the Hamiltonian takes the form of 

H = wl 2 . (65) 

The Kramers degeneracy of this Hamiltonian was discussed in Sec J2.51 Indeed, K symmetry here 
corresponds to = a x K. Thus by applying the following unitary transformation, 

e^ytev, v-±(\ 7), (66) 

reduces to the one in Sec J2.5( = K. At the same time, in this unitary transformation, q and 
H become 

q -)■ V^a z V = a y , H -»■ V*wl 2 V = wl 2 , (67) 

which are also the same as those in Sec l2.5l Thus the degeneracy here can be understood as a 
consequence of the generalized Kramers. 

4.2 e q = — 1 : pseudo-anti-hermiticity 

When e q = —1, the Hamiltonian becomes 

H = xa x + ya y + iza z . (68) 

From the general theorem 3 in Sec l3,21 H is expected to have a partner state with E and — E. 
Actually, the Hamiltonian has eigenvalues dry x 2 + y 2 — z 2 , thus the energy eigenvalues are paired. 
The eigenvalues are real when x 2 + y 2 — z 2 > 0, which corresponds to the model in quantum optics 
|28j. In this case, the constant energy surface in the parameter space is a one-leaf hyperboloid, 
H » . On the other hand, when x 2 + y 2 — z 2 < 0, the constant energy surface in the parameter 
space is described by the two- leaf hyperboloid, H 2 ' . Here we will consider the properties of the 
latter. 

The corresponding Hamiltonian takes the pure imaginary eigenvalues with opposite sign. 
Therefore, we transform H into iH , then deal with the following Hamiltonian, 

H = -xt x - yr y + zt z = Z %X V . (69) 

\ —ix + y -z J 

Here Tj (i = x,y,z) are "Pauli matrices" of SU(1,1): 

t x = if?x, T y = i<jy, t z = a z . (70) 

(See Appendix [Dj also.) As a consequence of the transformation H — > iH, the TR symmetry in 
the original Hamiltonian is converted into the particle/hole symmetry 

CHC- 1 =-H, C = a x - K, (71) 
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and the pseudo-anti-hermiticity becomes the pseudo-hermiticity, 

H^tjHt)- 1 , r] = a z . (72) 



Let us consider situations where all the eigenvalues of the Hamiltonian (|69p are real. Then, 
the parameters x,y, and z give coordinates on a two-leaf hyperboloid H 2,0 : 

x 2 + y 2 - z 2 = -r 2 < 0, (73) 

where r is a real positive constant. Since x and y are real, z is taken either z > r (upper leaf) or 
z < — r (lower leaf). The eigenvalues of the SU(1, 1) Hamiltonian ([69]) are given by 

E± = ±r. (74) 

On the upper leaf (z > r), the eigenvectors are given by 

|* + ) = -=L=( r + Z ), \<t,-) = ^=±=( " + ' X ), (75) 
y/2r(r + z) \ y-ix J ^/2r(r + z) \ r + z J 

which satisfy 

(<£±M0±> = (0±N^±>=±1- (76) 
Meanwhile, the eigenvectors of the hermite-conjugate Hamiltonian 

rrt I Z IX + y \ , . 

H ] = xt x + yr y + zt z = . , (77) 

\ ix — y —z J 

are 

1 / r + z \ .. 1 / —y — ix \ 

|^ + )) = ^=== , ¥>_ =-=== y , (78) 

A/2r(r + z) y -y + ix y y>2r(r + z) \ r + z J 

which satisfy 

{{<P±\V~ X \<P±)) = {{<P±Wz\<P±)) = ±1- (79) 
The bi-orthonormal bases, \4>±) and |y>±)), satisfy 

(^mbn)) = <W, ^ |0m)((^m| = |<M((<£+I + = I2, (80) 

and are related as 

-^±)=^|0 ± )=±|^ ± )). (84) 



;>= s=j J " y ~ fa )■ 1*'-)= * s ( 1. 



6 Similarly, for the lower- leaf (z < —r), the eigenvectors are 

1 ( -y - ^ \ , , 
y/2r{r - z) \ r-z J ' ~' s/2r(r - z) \ -y + ix 

and 

- vw^T) ( ' - ) ■ - ( ) ■ < s2 > 

They satisfy 

{(f>' ± \a z \cf>' ± ) = {{<p'±\o- z \cp'±)) = =Fl, {(t>'m\f' n )) = 5 mn . (83) 
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From the general arguments in Sec. 12.41 the particle/hole pair of \4> m ) is given by Crj 1 \(p m }}. 
With i] = a z and C = a x K, the particle/hole pair of \4>±) is given by 



Cr) 1 |y±)) = -ia y K\ip±) 



(85) 



Thus, the particle/hole pair of \4>+) is | </>_). Indeed, and \4>-) are eigenstates of H with 
opposite energies, EL = r and EL = — r. Though \(j>+) and \4>-) are not orthogonal in the usual 
sense, i.e. (0+|(/>_) 7^ 0, they are linearly independent. They are orthogonal in the sense of 
pseudo-inner product: 

(6-\<tM+) =0. 



Applying Crj 1 from the left to the both sides of (|84[) . f/l^i) = ±\tp±)), we have 

C|</» ± ) = ±Cr ? -V±))- (87) 
Therefore, Cr]~ 1 \(p±}} is equal to C\<fi±) up to sign. (See also (pOj) in Sec l2.4l ) Indeed, 

|0±) A |^>=C|^±). 

Thus, we find that the particle/hole pair of \4>±) is simply given by C\cp±) in the present case. 
The relations of bi-orthonormal bases, \4>±) and \(f±)), are summarized in Fig(TJ 



<l> + - 



9+ 



c 



Figure 1: The relations between the eigenvectors of the SU(1, 1) model. 



5 Random matrix class of 50(3, 2) model 

As a non-trivial realization of the generalized Kramers, we introduce 50(3,2) model. First 
consider K symmetry realized by 



k 

\ a x J 

which satisfies kk* = 1. Any arbitrary 4x4 matrix can be expanded as 

5 5 

a=l a<b=l 



(89) 



(90) 
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where I4 is 4 x 4 unit matrix, h, h a and h ab stand for complex parameters, y a denote 50(3,2) 
gamma matrices 

( a x \ ( <t v \ ( -ia-A / 1 2 \ fl 2 o\ 

(91) 

and 7^ are 50(3,2) generators constructed by 

7o6 = ^r[7a,76]- (92) 

The sixteen matrices, I4, -y a , y a b, amount to complete matrix bases that span arbitrary 4x4 
matrix. The imposition of K symmetry specifies h and h a as real parameters and h ab as pure 
imaginary parameters, 

h* = h, h a * = h a , h ab * = -h ab . (93) 
Thus, the Hamiltonian becomes 

5 5 
a=l a<b=l 

with real parameters, w, x a and x ab . We further impose Q symmetry 

which satisfies q* = —k~ l qk^~ 1 . According to two types of Q symmetry, e q = ±1, the Hamiltonian 
takes two different forms shown in Sec J5.ll and Sec J5.2t respectively. 

5.1 e q — +1: pseudo-hermiticity 

In this case, the Hamiltonian becomes 

5 

H = wU + Y J ^la. (96) 

a=l 

From the general theorem[2]in Sec ]3.2| this Hamiltonian is expected to exhibit Kramers degeneracy. 
The Hamiltonian is invariant under the time reversal symmetry 

QHG' 1 = H, (97) 

where 

cr x , 



@=l a * °} K, (98) 



with 9 2 = +1. 



This is the first nontrivial form of self-dual real split-quaternion matrix. See Appendix [C] for details. 
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The SO(3, 2) Hamiltonian (|96p is rewritten as 

\ x^i 2 + ix l n -x°i 2 J 

where xVj = xt x + yr,, + zt z . (Yli=i x%T i wm be abbreviated as xVj hereafter.) The eigenvalues 
of H are derived as E± = ± y/— (x 1 ) 2 — (x 2 ) 2 + (x 3 ) 2 + (x 4 ) 2 + (x 5 ) 2 . Note that we have two-fold 
degeneracy in the spectrum, which comes from the generalized Kramers theorem mentioned in 
SecEU 

Let us consider situations where all the eigenvalues of (|99h are real: 



{x ^f + (z 2 ) 2 _ ( x 3 ) 2 _ (x 4 )2 _ (x 5 )2 = _ r 2 < Q) (1Q0) 

where r is a real positive constant. (When x 4 = x 5 = 0, the SO(3, 2) model is reduced to two 
independent SU(1, 1) models.) The eigenvalues of the 50(3,2) model (f99|) are 

E± = ±r. (101) 

Here E + and £J_ are doubly degenerate, respectively. For x 5 > — r, the eigenvectors are given by 



" ^2r(r + x5) I (x 4 + zxV;)^ 



1 / (—x 4 + ixV, 

^(r + x 5 ) I (r + x 5 )^>, 



(103) 



where Q (a = ±) represent two-component spinors that account for double degeneracy. Take <p± 
as 

The hermitian conjugate of the Hamiltonian ([99p is 

x 5 lo x 4 lo — ixV„ 



Since 

f2 



\ X^l2 + IX Z T- —X°l2 J 

H ]2 = H 2 = r 2 l 4 , (105) 
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the eigenvalues of are also given by ±r, and the corresponding eigenvectors are|f| 

1 / (r + x 5 )4> a \ 



1 / (-x 4 + txVj)^ 



iV'ma) and |xmo)) are related as 



\lpma) f?|Vw) = a|Xma))- (HI) 

With (|1U2|) and (jllUp . it is straightforward to confirm that \ip ma ) and \xma)) indeed constitute 
the bi-orthonormal basis: 

(^ma|Xn/3» = (112) 
\^ma){{Xm a \ = \4>++) ((X++\ + |^+-> {(X+- 1 + |^-+) «X-+| + I ^— ) ((X— I = U- 

m,a=+,— 

With (jllip . the time-independent inner products induced by 77 are given as 

tymaMMnp) = {{Xma\v\Xnp)) = OL5 mn 5 a p. (113) 
Note the sign of these products depends on their "spin" directions. 



8 For x 5 < r, the eigenvectors are 

, _ 1 / (-x 4 + ix i T i ) 



IV'ia) and l^io) are related by the SU(1, 1) transformation 



<7± = — = =(~ £ 4 ± inn), (107) 



where g~ — g+ 1 . The eigenvectors of W are 

. , _ ( (-x 4 + ix'rl 

\X+al) — 



y / 2r(r - x 5 ) V -(r - 



^2r(r-a; 5 ) ^ -(sc + 
|x±«}) an d |x±Qf)) ar e related by the ST/(1, 1) transformation 



X- a ))= ,„/ J /">tu )• (108) 



Vr 2 ~ (z 5 ) 

where g'_ — <?+ _1 - 



= -== i_(-x 4 ± ixV/), (109) 
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As mentioned above, two-fold degeneracy is a consequence of the generalized Kramers theorem. 
The generalized Kramers pair of |Vw) is given by Qr]~ l \xma)) ■ Here, rj and are given by r\ = q 
(f95j) and (f98|) . Therefore, the generalized Kramers pair is derived as 



err 1 



X, 



v 
ia n . 



K\Xma)) = tt\A. 



(114) 



where a is the opposite spin of a; a = —a, i.e. (+) = — , (— ) = +. Thus, \tp ma ) and \ip m a) are 
indeed the generalized Kramers pair. Notice the "spins" a of the generalized Kramers pair are 
opposite to each other. Though they are not orthogonal in the ordinary sense, they are linearly 
independent. In the present case, they are orthogonal in the following inner product: 



i, 



0. 



(115) 



Since we have the relation (jllip which corresponds to (|24p with real eigenenergies, the generalized 
Kramers pair ®rj~ l \xma)) is equivalent to Q\ip ma ) (see Sec l2.3[) . Indeed, 



e 



Thus, we find the generalized Kramers pair is simply given by Q\ip ma ) 
The relations between \ip ma ) and \xmaf) are summarized in Figj2j 



¥±±- 



-X±i 



(116) 











X: 



Figure 2: The relations between the eigenvectors of the SO(3, 2) model. 

5.2 e q — — 1: pseudo-anti-hermiticity 

When e q = —1, the Hamiltonian takes the form of 

5 

H = iJ2 x0h ^- ( 117 ) 

a<b=l 

The Hamiltonian is invariant under the particle/hole symmetry 

CHC- 1 = —H, (118) 

where 

C (H9) 
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with C 2 = +1. From the theorem [3] in Sec ]3.2l we expect H has partner states whose energies 
are E and — E. Indeed, an explicit calculation shows the Hamiltonian has two paired states with 
energies (Ei, —E{) and (E 2 , -E 2 ) (in general E\ / .E^jE 

5.3 Realization as J = 3/2 SU(l, 1) quadrupole model 

In Ref.p], Avron et al. demonstrated that the 5 = 3/2 quadrupole Hamiltonian can be 
expressed by an SO (5) Hamiltonian. Here, we demonstrate how such arguments are generalized 
to the present non-hermitian case. The correspondences between SU(2) and SO (5); and SU(1, 1) 
and 50(3, 2) suggest that the SU(1, 1) spin 3/2 quadrupole Hamiltonian may be expressed by an 
50(3, 2) Hamiltonian. The SU(1, 1) spin Jj (i = x, y, z) are defined so as to satisfy 



[Ji, Jj] — iGijfcJ , 



(120) 



where denote the totally antisymmetric 3-rank tensor with e xyz = 1, and J 1 = (J x , J y , —J z ). 
With real 3x3 quadrupole coefficients , we introduce the SU(1, 1) quadrupole Hamiltonian as 



H(Q)= Yl Q i3J i J r 



(121) 



The SU (1, 1) quadrupole Hamiltonian is invariant under the SU(1, 1) spin flipping transformation 



Ji 



-Ji. 



(122) 



The five basis elements of Q %3 are taken as 
Qx 



i 



/ 1 


^10 

/ 1 

1 
^000 

which are orthonormal 




( i o \ 

0-10 




Qi 



(Qa,Qb) = -Tr(Q a Q b ) 



Jab- 



With the use of Q a , an arbitrary quadrupole matrix is expanded as 

5 



Q — ^ ^ X Q a , 



(123) 



(124) 



(125) 



a=l 



9 The expressions of E\ and Ei are rather lengthy, so we omit their explicit formulae 

(-1 

10 Note that Q5 is different from the traceless quadrupole matrix, i 



-1 I , used in Ref.[2]. 
2; 
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where x a are real, and the SU(1, 1) quadrupole Hamiltonian (|12ip is expressed as 

5 



where 



More explicitly, they are 



H(Q) = ^x a H(Q a ) 

0=1 

H{Qa)= (Q-T jj i J r 



(126) 
(127) 



H(Qi) = ^={J X , J z }, H(Q 2 ) = -L{J y , J z }, H(Q 3 ) = -L(jJ - J2), 

F(Q 4 ) = ^{ Jx, -V, ff(Q 5 ) = ^ + Jy + 2J 2 2 )- 
In particular, for J = 3/2, the SU(1, 1) spin is given by 4 x 4 matrices 



(128) 



/ V3i 



\ 



2i 
2i V3i 

^ o o v& o y 



-y/3 
-2 
\ 













( 3 











2 












1 











73 




' Jz ~ 2 








-1 





V3 

















-3 



The spin flipping operator 0' satisfying 



is given by 



e' 



a x 
a x 



K. 



\/2 



(129) 
(130) 

(131) 
(132) 



Since 0' satisfies 

e" = +i, 

0' can be regarded as the TR operator for 0' 2 = +1. (See Appendix [Dl for more details about 
the 577(1, 1) spin flipping operator.) Since H(Q a ) are the quadratic forms of the SU(1, 1) spins 
(|128p . they are invariant under the SU(1, 1) spin flipping operation. By substituting (|129p into 
(|128p . H(Q a ) are explicitly derived as 



H(Qa) = 7„, 



(133) 



with 7^ being 50(3,2) gamma matrices 



/ 

7i 



/ 

74 



ia x 








-ia x 





ih 


— H2 






/ 

72 



/ 

75 









-icr,, 



1 

7 3 



-1 2 
-1 2 



<x 2 
-a z 



(134) 
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Therefore, the J = 3/2 SU(1, 1) quadrupole model is expressed by 

5 

H = J2^'a, (135) 

a=l 

and the Hamiltonian is invariant under the TR transformation of 0' , [H, ©'] = 0. By rearranging 
the basis, 7^ (|134|) are transformed to the previous 50(3,2) gamma matrices j a (f£)Tj) . and 0' 
(|131|) is also to (|98p . Thus, we have shown that the J = 3/2 SU(1, 1) quadrupole Hamiltonian 
is equivalent to the 50(3, 2) Hamiltonian. 

6 Summary and discussion 

We explored the generalized Kramers degeneracy for 2 = +1 in pseudo-hermitian quantum 
mechanics. As the quaternions realize the TR operation for 2 = —1, the split-quaternions the 
TR operation for 2 = +1. We showed, by passing from the quaternions to split-quaternions, the 
following generalized theorems in pseudo-hermitian quantum mechanics: 

• If the system is invariant under the TR transformation 2 = +1 and also TR operator is 
anticommutative with the metric operator, the system has at least doubly degenerate states: 
the generalized Kramers pair. 

• When the system is invariant under the particle/hole transformation C 2 = +1 and also 
charge-conjugation operator C is anticommutative with the metric operator, the system has 
paired states with E and —E: the particle/hole pair. 

In both cases, the Hamiltonians necessarily possess the split-quaternion structure. We also identi- 
fied TR, particle/hole, and pseudo-(anti-)hermitian symmetries in the non-hermitian category pro- 
posed by Bernard and LeClair [27], and reconsidered the above theorems in view of non-hermitian 
random matrix. As a concrete example of the second theorem stated above, we investigated the 
SU(1, 1) model, and confirmed that the theorem indeed holds. Similarly, as an example of the first 
theorem, we introduced the 50(3,2) model. We confirmed that the 50(3,2) Hamiltonian is in- 
variant under TR transformation, and the TR symmetry brings double degeneracy to the 50(3, 2) 
model, exactly analogous to the Kramers degeneracy of the 50(5) model. The correspondences 
between the 50(3, 2) and J = 3/2 SU(1, 1) models are also clarified. 

As pointed out in Ref.[2], the structure of the original 50(5) model is related to instantons 
and twistor theory. Similarly, the present 50(3,2) model is related to split-instantons and 
twistor theory [30] . The present work was inspired by recent developments of the pseudo-hermitian 
quantum mechanics and the topological insulator in condensed matter physics. This work may 
hopefully be regarded as the first step of interplay between these two developments. It is intriguing 
to speculate realizations of the pseudo-hermitian quantum mechanics in condensed matter physics. 
We would like to pursue the issue in a future research. 
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A 577(1,1) and 50(3,2) 

Here we briefly review the non-compact groups, SU(1, 1) and SO(3, 2). 

The SU(1, 1) group consists of 2 x 2 matrices g satisfying the following relations, 

g ] a z g = a z , detg = 1. (136) 

Expanding g by its generators, T{ (i = x,y, z), g = 1 + i J2i 9i T i + ' ' " with real parameters Oi, we 
obtain 

r\a z - a z Ti = 0, trr* = 0. (137) 
Thus the generators of SU (1,1) are given by 



(138) 



where cij (i = x,y,z) are the standard Pauli matrices for SU(2) group. The SU (1,1) Pauli 
matrices, n = (t x ,t v ,t z ) = (icr x ,ia y ,a z ), satisfy the following relations: 

VxTiCTx' 1 = t* , (139a) 
OyTiOy' 1 = -Tj*, (139b) 
oyrjCTz -1 = tJ. (139c) 

The SO (3, 2) group is linear transformations with unit determinant acting on a five dimensional 
vector (x\, X2, x$, X4, X5) and preserving the following norm, 

-xj-xl + xj + xl + xl. (140) 

Its element is given by a 5 x 5 matrix G which satisfies, 

^2 G acG b dVcd = Vab, detG = 1, (141) 

cd 

with rjab = diag(— 1, —1, 1, 1, 1). Expanding G by its generators M ab (a, b = 1, 2, • • • , 5), we find 

(M ab ) ceVed + (M ab ) deVce = 0, (M ab ) cc = 0, (142) 

where we use the convention in which the repeating indices are summed. These relations are met 
by the following M ab , 

(M ab ) cd = i(5 ca r] bd - 5 cbVad ), (143) 
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which satisfies 

[M ab , M cd ] = -i(ri ac M bd - Vbc M ad - Vad M bc - Vbd M ac ). (144) 

Now consider the spinor representation of 50(3, 2). The spinor representation is given by the 
4x4 matrices 7 a with anticommutation relations 

{7a, 7b} = 2r?afe- (145) 
The anticommutation relation is obtained by 

( a x \ ( a y \ ( -ia z \ 

and the 50(3, 2) algebra (|144|) is realized by 

M ab = lab ^^ci a , lb \- (147) 

The gamma matrices (|146p and generators ()147|) of 50(3, 2) satisfy the following relations, 

k^ak" 1 = 7a*, k^ahk" 1 = -j ab * , (148a) 

C7aC~ 1 = 7a*, C7 ab c _1 = -7ao*, (148b) 

97ag _1 = 7a f , qiabQ' 1 = lab 1 , (148c) 



where 



a x a, a z 

J' I a y j ' I <r 2 



(149) 



B Quaternion and split quaternion 

The quaternion (1, ei,e2, 63) are defined so as to satisfy [21 



(ei) 2 = -l, (e 2 ) 2 = -l, (e 3 ) 2 = -l, 

e^ej = -e^e* (i / j), eie 2 e 3 = -1. (150) 

The "imaginary" quaternions are realized as Pauli matrices as 

(ei,e 2 ,e 3 ) = (-ia x ,-ia y ,-ia z ). (151) 

The split-quaternion algebra (1, q\, q2, 53) is simply obtained by flipping two signs of squares of 
quaternions: 

(qi? = +1, te) 2 = +1, te) 2 = -1, 

QiQj = -QjQi {i / j), 919293 = -1- (152) 
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The split-quaternions are realized by non-hermitian matrices as 

Q2, 93) = (iT x ,iT y ,iT z ) = ( c x , -a y , ia z ), (153) 

where T{ (i = x,y,z) denote the SU(1, 1) "Pauli matrices" (|138[) . They satisfy 

[Ti,Tj] = 2ie ijk T k , {Ti,Tj} = -2rjij, (154) 

where is the three rank antisymmetric tensor with e xyz = 1, while rjij = diag(+l, +1, —1) and 
t 1 = (T x ,T y , —t z ). Replacing the imaginary unit i in a y (|15ip with three imaginary quaternions, 
the Pauli matrices are "enhanced" to yield 5*0(5) gamma matrices: 




(155) 



It is straightforward to see that (| 155j) satisfy {^ a , jb} = 25 a b. By applying such substitution in the 
case of split-quaternions, we obtain 4x4 non-hermitian gamma matrices of 50(3,2) (|146p . The 
correspondence can also be naturally understood by noticing isomorphism of groups: SU(2) ~ 
USp(2) and 50(5) ~ USp(4); SU (1,1) ~ Sp(2,R) and 50(3,2) ~ Sp(4,R). 



C Definitions and relations for split-quaternions 

We introduce the terminology for split-quaternions in the same spirit for quaternions (see 
Refs.[TJ[23] for instance). The split-quaternion generally takes the form of 

q = cq + ctqi, (156) 

where c and q (i = 1,2,3) are complex numbers. There are three types of conjugation for split- 
quaternion: The complex conjugate, split-quaternionic conjugate, and split-quaternionic hermitian 
conjugate, which are respectively defined by 

q* = c*q + c*qi, (157a) 
q = cq -Ciqi, (157b) 
q t = {q*) = c*qo ~ c* Q i- (157c) 

Such conjugations have the following properties: 

(<?i • Q2)* = qi* ■ q2*, (qi ■ 92) = <fe ■ qi, (li ■ <?2)* = Q2 t ■ q^- (158) 
With the matrix realization (|153|) . split-quaternion (|156p is expressed as 

. . / c + ic 3 -c x +ic 2 \ , K ■. 

q = c + CiiTi = c- c\a x - c 2 cr y + c 3 ia z = . . (159) 

\ -ci - ic 2 c- ^C3 / 
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Correspondingly, the three kinds of conjugate (|157p are 

q = c + Ci iTi = c - c x a x - c 2 a y + c 3 ia z = , . „ „ , (160a) 



-1 



g = c - CjZTj = c + ciCTz + c 2 a v - c 3 ia z = . , (160b) 

\ ci + ic 2 c + zc 3 / 

ql = c - Ci in = c + Cl a x + Gjffy - Cgtff, = : . (160c) 

Due to the non-hermitian property of the split-quaternions, t£ = (—T x , —T y , T z ), the split-quaternionic 
hermitian conjugate (|160cp does not coincide with the ordinary definition of the hermitian conju- 
gate 

g t = c * - c * iTi \ = c *- C la x - c* 2 a y - c%ia z = ( °* ~ ^ ~f + ) . (161) 

\ -C* - IC* 2 C* + icl J 

(In the quaternion case, the quaternionic hermitian conjugate coincides with the ordinary hermi- 
tian conjugate.) 

The real split-quaternion is defined as 

r , i ( w + ix 3 -x 1 +ix 2 \ 

q = wq + x qi = , . „ . „ , (162) 

\ —X — IX w — IX J 

where w and x % (i = 1,2,3) are real numbers. The necessary and sufficient condition for the real 
split-quaternion is given by 

qi = q. (163) 

An M x M split-quaternion matrix (2M x 2M matrix in the usual sense) Q is defined as a matrix 
whose matrix elements are split-quaternions: 

(Q)lJ = QIJ, (164) 

where I,, J = 1,2, ■■■ , M. The complex conjugation, split-quaternionic conjugation, and split- 
quaternionic hermitian conjugation of Q, are respectively defined as 

(Q*)lJ = qiJ*, (165a) 
(Q)lJ = qJI, (165b) 
{Q X )lJ = qjl X - (165c) 

We call Q the "dual" of Q. The split-quaternionic hermitian matrix is a split-quaternion matrix 
that satisfies 

Q X = Q. (166) 
Unlike quaternion matrix, the split-quaternionic hermitian matrix is not hermitian in the usual 

/ hj — — ZX^ — X^ \ 

sense. For instance, q = w + ix l qi = w — x l Ti = ( 1 „ , (with real numbers 

\ — ix + x w + x / 
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w , x 1 , x 2 , x 3 ) is split-quaternionic hermitian, but not hermitian in the usual sense. A real split- 
quaternionic matrix refers to the matrix whose components are real split-quaternions qjj, 

(Qlu = qh, (167) 

and then Q r satisfies the relation 

Q t = Q. (168) 

The self-dual real split-quaternion matrix is defined as a split-quaternion matrix that satisfies 
both (fTMl) and (fTBHjl : 

Q = Q t = Q. (169) 

Thus, the split-quaternionic hermitian real split-quaternion matrix is equivalent to the self-dual 
real split-quaternion matrix. (The condition Q = Q is the self-dual condition.) The terminology 
"split-quaternionic hermitian real split-quaternion" is rather clumsy, so we use "self-dual real 
split-quaternion" instead. Such self-dual real split-quaternion matrix generally accommodates 
the generalized Kramers degeneracy for 9 2 = +1 In low dimensions, the self-dual real split- 
quaternion matrices are given by 



M = 1 



M 



Qi 




W + X 



X + IX 





w + x 5 

1 ■ 9 

-X + IX 



-X 



IX 



IX 



■ S 1 • 2 \ 

x" — IX x — IX \ 

x 1 + ix 2 x 4 + ix 
w — x 5 




w — X 



(170) 



where uu,x , • • • , x 5 are real parameters. With the 50(3, 2) gamma matrices (|9ip . Q2 is concisely 
represented as 



w 



1 4 + j^X°7a- 



(171) 



a=l 



Qi and Q2 are exactly equal to the matrices (|65|) and ([96]) . respectively. They have both K and 
Q symmetries {e q = +1), and their eigenvalues are 



M 
M 



1 : 

2 : 



w ± yj-{x 1 ) 2 - (x 2 ) 2 + (x 3 ) 2 + (x 4 ) 2 + (x 5 ) 2 , 
with double degeneracy (of the generalized Kramers). 



(172) 



D Spin nipping operators and quadrupole Hamiltonians for low 
SU (1, 1) spins 

The SU(1, 1) algebra is given by 

•/,••// /(„/,./ A '. (173) 

Real split-quaternion matrix does not accommodate the generalized Kramers degeneracy in general. The 
split-quaternionic hermitian condition has to be imposed as well. 
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where denotes a totally antisymmetric tensor with e xyz = 1, and J 1 = (J x , J y , —J z ). Explicitly. 

\Jxi Jy\ — iJzi [Jyi Jgj — iJx) \Jzt Jx\ — iJy (174) 

From SU(2) spins S x ,S y ,S z , the SU (1,1) spins are constructed with the identification 

J x = iS x , Jy = iSy, J z = S z . (175) 

Note that J x is pure imaginary; J y and J z are real. The magnitude of SU(1, 1) spin J is defined 

as 



J 2 X -J 2 V +J 2 Z =J(J + 1). 



For instance, 



_ 1 _ 1 _ 1 

J — 1/2 : J x — 2 Tx ' ~ 2 Tyi ^ z ~ ~2 Tz ' 



(176) 
(177a) 



/ i 
i i 
^ i 



j i \ , Jy = —f=. 



1 

V2 



/ o 


1 0^ 




( 


1 


\ 


-1 


1 












V 


-i oy 




{ 





-l) 



The SU(1, 1) spin flipping operator H is defined so as to satisfy 

J% ^ E*J%£* — J%. 

Express S = ?7 ■ K, and U satisfies 

UJ X U~ 1 = J X , UJyU- 1 =-Jy, \JJ Z \J- X = -J Z . 

Here, U is a unitary matrix given by 

U = (-i) 2J e nJ *, 

where the factor (— i) 2J is added for convenience. Consequently, H is given by 

S = (-i) 2J e wJ * ■ K, 

which satisfies 

r? = +1, 

13 



independent of the magnitude of the SU(1, 1) spin I. For low SU(1, 1) spins, 

J =1/2 : E 1/2 = a x -K, 

( 1 

J = 1 : Si = 



V 



1 I K, 

1 



J - 3/2 : H 3 / 2 



a x 
a x 



• K. 



(177b) 

(178) 

(179) 

(180) 

(181) 
(182) 

(183a) 
(183b) 

(183c) 



12 As discussed in Sec l2.6l the square of the SU(2) spin flipping operator takes —1 for half-integer spins, while 
the 517(1, 1) spin flipping operator yields +1 even for half-integer 5(7(1, 1) spins. 
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Note H1/2 is equal to the charge conjugation operator C (fTTj) of the 577(1,1) model, and H 3 / 2 
the time-reversal operator 6' f 1 1 3 1 [) of the 50(3, 2) model. Similarly, for low 577(1,1) spins, the 
quadrupole Hamiltonians (|127p introduced in Sec J5.3l are given by 



J =1/2: H(Q a )=0, 



J=l: fl-(Qi) 



(184a) 



H(Q 3 ) 




( 1 \ 



#(Q 2 ) 



77(Q 4 ) 




/ 1 

0-2 
1 
(184b) 



E V, T, and (J operators for SU (1,1) and 50(3,2) models 



For non-hermitian Hamiltonians, there is a systematic procedure for constructing a metric 
operator 770, and symmetry operators which commute with the Hamiltonian |31| I32j . In this 
Appendix, we briefly review the procedure and apply it to the ££7(1,1) and 50(3,2) models, 
respectively. 

For a non-hermitian Hamiltonian H, the Schrodinger equation is 



H\<j> n ) = E n 



rf\y n ))=E* n W, 



(185) 



where the eigenvalues E n are complex in general. As discussed in SecEjJ the eigenvectors, \4> n ) 
and \<p n )), give a bi-orthonormal basis [23|, I22j satisfying the orthonormal and complete relations 



((<Pm\<j>n) = (4>m\<Pn)) 



\<t>n){(<Pn\ = ^ l^n))(0n| = 1. 



The Hamiltonian and its hermitian conjugate are expanded as 

n n 

A non-hermitian Hamiltonian is called pseudo-hermite when it satisfies 

ift = rjHri' 1 , 

where r\ is hermitian and called a metric operator. A state given by 

\4>n)' = rf x \(p n )), 



(186) 



(187) 



(188) 



(189) 



is an eigenvector of H with an eigenvalue E* as seen from (|185p and (|188|) . Thus the eigenvalues 
of the pseudo-hermitian Hamiltonian are classified into two types. One is a set of real eigenvalues 
and the other is a set of complex conjugate pairs. (For more details, see [T7].) 
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Suppose all the eigenvalues of a non-hermitian Hamiltonian are real, i.e. E* = E n . Following 
Ref.[32], a metric operator 770 is constructed as 



N N 



*7o 1 = X>>^l' 7 » = Z)lw))((^l- (190) 
1=1 1=1 

Here note that the metric operator satisfying (|188|) is not unique. t]q in the above generally 
depends on parameters of the Hamiltonian, while we may have a constant metric operator for 
some models. For example, see Sec lE.ll and Sec lE.2l 

Let us now define the following operators V, T, and C, 

P\tp n )) = (-l) n+1 |<M, (191a) 
T\4>n) = \<Pn)), (191b) 
C^n) = (-l)^ 1 !^). (191c) 

These operators were originally introduced in analogy with parity, TR, and charge conjugation, 
respectively, however they are not, in fact, directly related. With the orthonormal and complete 
conditions (j!86|) . they are explicitly written as 

N 



^ = E(- 1 ) /+1 i^)^i' ( 192a ) 
1=1 

N 

T=Y,Wl))Kiw\, (192b) 
1=1 

N 

C = E(- 1 ) m |^)((^| ) (192c) 



which are respectively hermitian, anti-hermitian 13 !. and pseudo-hermitian 

V ] = V, (193a) 
T f = T, (193b) 
C^tioCtiq 1 . (193c) 

From (|192p . VT and CVT operators are given by 

N 



PT = J2(- i y +1 \h)K((vil (194a) 
1=1 

N 

CVT = Y,\&)K((w\, (194b) 



1=1 



13 The terminology "anti-hermitian" usually refers to operator whose hermitian conjugate is equal to the minus 
of the original, but here, the terminology "anti-hermitian" refers to hermitian operator that anticommutes with 
imaginary unit. 
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and 



They are pseudo-antiunitary 



V T\<p n ) = (-l) n+l \ ( f> n ), (195a) 
CVT\(f>n) = \4>n)- (195b) 



{VT? = Vo (PT)- 1 Vo 1 , (196a) 
(CVT? = V o(CVTr 1 ^ 1 . (196b) 



It is readily seen that these operators satisfy 



[H, VT] = 0, 

[H,C] = [H,CVT]=0, 

[VT,C]=0, 

[H,V}^0, (197) 



and 



(VT) 2 = C 2 = (CVTf = 1, 

V 2 / 1, 7~V 1- (198) 

As realized in the first and the second lines of (|197|) . the Hamiltonian always displays "VT 
symmetry" and "C symmetry " (or ll CVT symmetry ") with respect to the VT and C operators 
constructed above. 

E.l 517(1,1) model 

With the bi-orthonormal bases (|75p and jTHJ), we construct the metric operator ryo, and V, T, 
and C operators for the SU(1, 1) model. Prom (|190p . 

r j0 = \< P+ ))(( V+ \ + \^))((^\ = -(, Z - lx -y), (199) 

r \ix — y z I 

from (fTMj) . 

V = \ ( f> + )(<P + \-\<P-)(^\=a z , 

T (\in \\//rn* I -L \,n h K 1 ( {v + zf + {%X + V f -2{r + z) V \ 

T= (\<P + )) (M + \<P-)) K= WT - ) \^ _ 2{r + z)y {r + z? + {lx + y? )- K > 

C = \<P + )(^ + \-\^)((^\ = -( .\ - lX - y ) = 1 -H, (200) 

r \ —ix + y —z I r 
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with H ([Ml), and from (fTMjl . 

pr = (|^;|H*->«<l)■^• = ;^7 V^^ (r + ^t < t + ! ' ,2 , A* 

2r[r + z) \ 2(r + z)y — (r + z)—(ix — y) I 



1 / (r + z) 2 - (ix + y) 2 2i(r + z)x 

2r(r + z) \ —2i(r + z)x (r + z) 2 — (ix — y)' z 

(201) 



The metric operator (|199p is different from 77 = cr 2 used in SecJH This "discrepancy" stems from 
the non-uniqueness of the metric operator. 

E.2 50(3,2) model 

For 50(3,2) model, with use of the bi-orthonormal basis |Vw) (|102p and \xmaf) flUOp . the 
metric operator is constructed as 

% = dX+a»((X+a| + IX-o»((X-o|) 

a=+- 

1 / (r + x 5 ) 2 — (x 4 — ixV/)(x 4 + ix^Tj) — i(r + x 5 )x*(rj — r/) \ 

2r(r + x 5 ) I — i(r + x 5 )x*(rj — r|) (r + x 5 ) 2 + (x 4 + ix l rj)(x 4 — ixHj) J 

(202) 

Similar to the case of SU (1, 1) model, the metric operator (|202p is different from the one given by 
([MD . ([95]) is anticommutative with 6 (JSHJ), while (p02|) is commutative with 9. From (fT92|) . V, T 
and C are derived as 

V= ]T Cl^+aXV'+al-l^-aX^-al) 
a=+- 

1 / (r + x 5 ) 2 -(x 4 -ixVi)(x 4 + ^rj) (r + x 5 )(2x 4 - ix^n + t})) \ 

~ 2r(r + x 5 ) y (r + x 5 )(2x 4 + is* (75 + r\)) -(r + x 5 ) 2 + (x 4 + ixV;)(x 4 - ia^Vj) J ' 

T= Y dX+«))((X+al + IX-a»((X-oJ) • K 



Q = 



1 / (r + x 5 ) 2 + (x 4 - ixV/)(x 4 - ix J r*) i(r + x 5 )x l (r* + tJ) 

2r(r + x 5 )l i(r + x 5 )x*(t* + r/) (r + x 5 ) 2 + (x 4 + zxV/)(x 4 + ix j r*) 



■K, 



C= Y, (|^+a)((X+a|-|^-a)((X-a|) 

-H, (203) 



1 / X 5 l2 X 4 lo — ix l Tj\ 1 



L2 

,4 _l 



r \x -\-ix l Ti — x s i2 / r 
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with H ([M]), and VT and CVT are 

PT= £ (l^+a)((X^|-|^-«>«X-«l)-^ 

1 / (r + x 5 ) 2 — (x 4 — ix l Ti)(x A — ix^T*) (r + x 5 )(2x 4 — ix J (Tj — r*)) 



2r(r + x 5 ) I (r + x 5 )(2x 4 + ix i (r i -r*)) -(r + x 5 ) 2 + (x 4 + ixV;)(x 4 + ixH* 

CVT= £ (l^+a>«X^I + l^-a><(X-al)-^ 
a=+ - 

1 / (r + x 5 ) 2 + (x 4 - ixVj)(x 4 - ix j Tj) i(r + x 5 )x J (rj + r*)) 



2r(r + x 5 )l i(r + x 5 )x J (rj + r*) (r + x 5 ) 2 + (x 4 + ixVj)(x 4 + ix J T*) y 

(204) 

F Level crossing point and monopole structure in non- hermit ian 
Hamiltonians 

F.l S77(l,l) model and U(l) monopole 

Singularity of phase of eigenstate generally reflects the non-trivial topology in phase space 
[3111]. For instance, the crossing point of two energy levels of the SU(2) model is called a diabolic 
point (an isolated point), which brings U(l) holonomy in phase space [33] . In the SU(1, 1) model, 
the eigen-energies are given by E± = ±r with r = \J z 2 — x 2 — y 2 , and the level crossing point 
E + = E_ is achieved when r = 0. In the SU(2) model, because of the Euclidean signature, the 
condition, r = \J x 2 + y 2 + z 2 = 0, is met only at the point x = y = z = 0. Meanwhile, in the 
SU(1, 1) model, the signature is hyperbolic, and the condition is satisfied on the surface 

x 2 + y 2 - z 2 = 0. (205) 

In such a case, r = point is called the exceptional point. Around the exceptional point (r ~ 0), 
the upper and lower energy eigenvectors {z > r) (|75p behave as 



1 / z \ . . 1 I y + ix 



V2rz \ y — ix J \/2rz \ z 

and 



(206) 



V2rz \ -y + ix / V2rz \ z 



They are degenerate, as found 



~ e **|0 + >, 

\<P-)) ~ ~e lx W + )), (208) 



where tanx = -• Then, the normalization condition is not satisfied but 



v 



(0 + |^ + )) = (0_|^_))~f ~0, (0+|p_)) = O, (0_|p+))=O. (209) 
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The holonomy of the exceptional points is 

A + = -i((<p + \dct> + ) = -i(<t>+\d<p + )) = -i(<f> + \a z \d(t> + ) = dx l Af, (210) 

where 

A t = o ( "T ^ iJ^, (211) 
2r(r + z) 

which is the U(l) monopole gauge field in a hyperbolic space [34} 135] . The corresponding field 
strength is calculated as 

= diAf - djAf = -^e iik x k . (212) 
Similarly, the holonomy for the negative energy state is evaluated as 

A' = -%({tp_\d4J) = -i(cj)_\d(pj)) = -i{cj)_\{-a z )\dcj)J) = dx*A-, (213) 
where the insertion of —a z is in accordance with the normalization (|76h . and 



The corresponding gauge field strength is 

F~ 3 = d t A~ - d 3 A~ = ~^e ljk x k = -F+. (215) 

The field strength diverges at the exceptional point. 

On the lower leaf {z < — r), the holonomies are derived as 

A' + = -iy + \d<j>' + ) = -i(cf>' + \d<p' + )) = -i{<f>' + \(-a x )\d^' + ) = dx l A'+, (216a) 
A'_ = -i((ip' _\d(j)'_) = -t(0'_|d^'_» = -i((f)' -\a z \d(/)'_) = dx'A'f, (216b) 

where 

A[ + = -AT = -—^—e^xK (217) 
2r(r — z) 

As found above, the holonomy of the exceptional point in the SU(1, 1) model is regarded as 
the gauge field of hyperbolic U(l) monopole, and the monopole charges for the upper and lower 
energy states are opposite. Such effect is similar to the U(l) monopole holonomy of the diabolic 
point in the SU{2) model [33]. 



F.2 50(3,2) model and 577(1, 1) monopole 

Degeneracies in energy levels generally bring non-abelian holonomy in the parameter space 
|36j. For instance, the 5*0(5) Hamiltonian (Luttinger Hamiltonian |37j ) has SU(2) holonomy 
which is crucial for the spin-Hall effect |38| . Here, we consider what kind of holonomy could 
emerge in the 50(3,2) model. The energy levels of the 50(3, 2) model are ±r, and the level 
crossing point is at r = 0, namely 

( x l)2 + (x 2 )2 _ (x 3 )2 _ (x 4 )2 _ (a . 5) 2 = , (21g) 
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Near the exceptional point (r ~ 0), the upper energy and lower energy eigenvectors behave as 



and 

I w 1 I x5 a i I w 1 ( -(x 4 - ix l T})4>a \ , , 

V2rx 5 y [x + ix l Tl)(p a J V2rx 5 \ x°<p a J 

Then, at the exceptional point, \ip+ a ) and \ip-a) are related by the SU(1, 1) gauge transformation 

— x — 1 $ Ti , and \x+a)) and |x— a}) are by — - — r ' , and the normalization conditions are no longer 
satisfied: 

<<M*±a'» ~° ^5 ~ 0, ftMX T a'» = 0. (221) 

For x 5 > — r, the upper energy degenerate eigenvectors bring the following holonomy: 

A + = -i((x+|dV+> = -^+|dx+» = dsVW^ (222) 

where 

= - 2r(r + g5) ^^ ^ = °- (223) 

Here, fx, v = 1,2,3,4, and rj^i are the "split"-'t Hooft symbol defined by rj^ui = e^ u i + rj^rj^ — 
VviVv'i with riftv = diag(+, +,—,—), and e^ u i is a totally antisymmetric tensor, ei23 = 1 and 
^fivi = for /x = 4 or z/ = 4. The corresponding field strength 

^ = ^+ - + i[A+,A+], (224) 

is derived as 

1 - 



*3 = --PS = ^(r + x 5 )A+ (225) 
The gauge field strength has the singularity at the exceptional point. Similarly, A~ is given by 

A- = -t«x-|#-> = -i(^_|dx-» = cte>W-</>, (226) 

where 

^ = ~ 2r(r + ^) r? ^ V ' ^ = °' (227) 
with 77^ = e^j - + ^i^4- The corresponding field strength is 

^5 = = + (228) 
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With a different gauge choice (x 5 < r), the holonomy is calculated as 

A' + = -i(( X '+W + ) = -i(/ + |dx+» = dx a <fta,A't<f>, (229) 

where 

< = - 2r(r _^ 5) i^ < = 0- (230) 
The corresponding field strength is 

r iiv ~ r 2 x ^ ^ r 2 " /•» 2r 2 ^ ' 

F'+ =-F'+ =-^(r-x 5 )<- (231) 
Similarly, the lower energy degenerate states bring the holonomy 

A'~ = = = dx a ^a z A'~^ (232) 

where 

A » = ~2r(r - x^)^^ T% ' A ' 5 = °' ( ^ 233 ' ) 
The corresponding field strength is 

F'^ = -F^ = ±(r-x 5 )A'+. (234) 

The SU(1, 1) gauge transformation relates A + and ^4 /+ as 

o z A'+dx a = ^(a 2 A+(ix a ) 5+ - ig j + a z dg + , (235) 

and their field strengths as 

<T z F , ab = g t + (<T z F ab )g + , (236) 
where g± are given by (|107p . Similarly, A~ and A' - are related by the SU(1, 1) transformation, 

a z A'~dx a = gl{a z A-dx a )g^ - igla z dg_, (237) 

and 

* z F'- b = gl{a z F- b )g_. (238) 

Thus, the exceptional points of the 50(3, 2) model act as the SU(1, 1) monopole with opposite 
charges for the upper and lower energy states. Such non-compact gauge group monopoles have 
been introduced in the context of the non-compact Hopf maps [35] . 
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